Strong-Segregation Theory of Bicontinuous Phases in Block Copolymers 
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We compute phase diagrams for A n B m starblock copolymers in the strong-segregation regime 
as a function of volume fraction <j>, including bicontinuous phases related to minimal surfaces (G, 
D, and P surfaces) as candidate structures. We present the details of a general method to compute 
free energies in the strong segregation limit, and demonstrate that the gyroid G phase is the most 
nearly stable among the bicontinuous phases considered. We explore some effects of conformational 
asymmetry on the topology of the phase diagram. 



I. INTRODUCTION 



Block copolymers (BCPs), comprising chemically dis- 
tinct polymers permanently linked together, are inter- 
esting because of the diverse array of ordered phases to 
which both polymer theory and experiment have been 
directed. The phase behavior of diblock copolymer 
melts is a competition between the entropic tendency to 
mix the two species into an isotropic melt, and an en- 
ergetic penalty for having unlike species adjacent, which 
induces transitions into ordered phases of many symme- 
tries, depending on the topology and composition of the 
polymers. Near the order-disorder transition (weak in- 
compatibility) entropy effects dominate, and the indi- 
vidual polymers retain (within mean field) their Gaus- 
sian coil conformation through the transition, while 
at much higher incompatibilities the chains are strongly 
stretched. It is this strongly stretched regime which we 
address here. 

Leibler developed the first complete theory of ordered 
phases in BCP melts , and predicted the by-now clas- 
sical phases of lamellar (L), cylindrical (C) and spherical 
(S) symmetry using the random phase approximation to 
derive an effective Landau free energy in terms of com- 
position modulations in Fourier space. The strong segre- 
gation regime was studied by Helfand and co-workers || 
and Semenov || , who predicted the same series of phases 
with increasing asymmetry, denoted by the fraction (f> of 
polymer A in an A—B diblock. (In this work we always 
use A to denote the minority block) . This treatment bal- 
ances the stretching energy of a polymer chain with the 
interfacial energy between A and B regions. By assuming 
an incompressible melt, minimization of the free energy 
gives a preferred domain size which scales as TV 2 / 3 , where 
N is the degree of polymerization. 

In the strong segregation limit the free energies of all 
microphases scale the same way with chain length and 
interfacial tension, so the phase boundaries become in- 
dependent of the strength of the repulsion \ between A 
and B monomers and depend only on the composition /. 
Semenov's calculation in effect gave a lower bound to the 



free energy of the L, C, and S phases because the phases 
he constructed did not fill space, but were micelles of the 
corresponding topology fl. This approximation treats 
the A — B interface and outer block surface as having 
the same circular or spherical shape, and is sufficient for 
understanding the qualitative aspects of the transitions 
between the phases. 

Experiments followed the theories of Leibler and Se- 
menov and quickly discovered a new phase, |g|-|l0[|, orig- 
inally thought to be ordered bicontinuous double dia- 
mond (here denoted D), of Pn3m symmetry, but recently 
shown to be of Ia3d symmetry [pl]-p^| and related to 
the minimal surface known as the gyroid (G). |l4j The G 
phase occurs for compositions between those of the L and 
C phases, can occur directly from the disordered phase 
upon increasing the incompatibility ^iV, and is found to 
be unstable to the L or C phases at high enough 

Although several groups attempted to describe this 
transition theoretically, |l5|-[l7| using variations on 
Leibler's theory, the first successful theory is due to Mat- 
sen and Schick [|l8|, who developed a method for com- 
puting the free energy of any crystalline structure by ex- 
panding the partition function in the basis functions for 
the symmetry of the desired mesophase, rather than the 
Fourier mode expansion of Leibler. They found a stable 
gyroid phase for 11.14 < %A < 60, where the upper limit 
was determined by extrapolation from the phase bound- 
aries at lower %A. |l9| This was followed by careful ap- 
plication of Leibler's method, pO| , ^l| to include higher 
harmonics and calculate the stability of the G phase in 
weak segregation analytically. 

Roughly concurrent to the calculations of Matsen and 
Schick, methods were developed to calculate the free 
energy of essentially arbitrary structures in the strong 
segregation regime (\N — ► oo). 0,^2). These meth- 
ods use the results for polymer brushes, |q]23|], supple- 
mented by an ansatz about the geometry of the rele- 
vant phase and an assumption about the chain paths. 
Olmsted and Milner assumed straight paths through the 
A — B interface and locally specified the volume fraction 
per molecule, |?],|24|,^5|] , while Likhtman and Semenov re- 
laxed the assumption of straight paths 123] but enforced 
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the constraint of constant <j> per molecule only globally. 
The former approach corresponds to an upper bound on 
the free energy (see below), while it is not clear that 
the Likhtman-Semenov calculations corresponds to any 
bound, or indeed to any systematic approximation, be- 
cause the local constraint of constant composition is re- 
laxed. By comparing upper bounds between bicontinu- 
ous, C, and L phases (obtained for the cylindrical phase 
by assuming hexagonal symmetry and imposing straight 
paths), we showed that the bicontinuous phases are un- 
stable, when comparing upper bounds, to the L and C 
phases. Later, Xi and Milner extended this work to cal- 
culations with kinked polymer paths, and found an upper 
bound to the hexagonal phase which lies very close to the 
lower bound using round unit cells. |f2(| 

Experiments have found an additional phase at x va l~ 
ues between the G and L phases ps|] , a hexagonally- 
perforated lamellae (HPL) phase, which consists of ma- 
jority lamellae connected through a minority matrix by 
hexagonal arrays of tubes. |29| The stacking has been 
suggested to be ABC ABC |J or ABAB (28). Theo- 
retical attempts to justify this phase have failed in both 
the strong segregation limit, where Fredrickson chose a 
catenoid as a candidate base surface; (30| and in the 
weak-segregation limit by self-consistent field calcula- 
tions 0]. Recent experiments have shown that the 
HPL phase is not an equilibrium phase in diblock melts, 
but may be metastable. 

Here we present the calculations of Rcf. |Q] in more 
detail. We show that the G geometry is the most stable 
of the candidate bicontinuous phases, followed by the D 
and P geometries, and that the G phase can be stable for 
block-copolymers with sufficient conformational asym- 
metry. The outline of this paper is as follows. In Sec- 
tion || we present the formalism for calculating the free 
energy in general geometries. In Section [1II| we present 
the results for the classical diblock topologies (lamellae, 
cylinders, and spheres), extended to include non-round 
unit cells and, in the case of the cylindrical topology, 
kinked paths. In Section |v] we present the free energy 
for a generic "saddle" wedge, which is representative of a 
generic bicontinuous structure as a pie-shaped wedge is 
representative of cylindrical phases regardless of packing. 
We then introduce the geometry necessary for calculat- 
ing the free energy of the P, D, and G topologies. In 
Section ^ we present our results for both symmetric and 
non-symmetric stars, and we conclude in Section VI. 



II. 



STRONG SEGREGATION THEORY IN 
GENERAL GEOMETRIES 



A. A single wedge 

We first recall some results for polymer brushes in 
strong segregation under melt conditions, and then show 
how to apply this to a general geometry. We consider a 



melt of star A nA —B nB copolymers , comprising n A arms 
of yl-blocks of mean square end-to-end distance R A and 
similarly for the 5-arms. The volume fraction <\> of A 
material is 1331 



\ ^ A 



(2.1) 



where Q = n A Q A + n B fl B is the total chain volume, and 
f2 A and fl B are the volumes of single A and B arms. 
Our calculations are appropriate for strongly-segregated 
chains, for which interfaces are sharp on the scale of mi- 
crophase lattice constants. In strong segregation the free 
energies of all microphases scale the same way with chain 
length and interfacial tension, so the phase boundaries 
become independent of the strength of the repulsion \ 
between A and B monomers. 




R 



FIG. 1. A typical wedge filled with an A — B star copoly- 
mer, with n A — 2 A arms and n B — 1 B arm. 

Consider an elementary wedge, as in Figure [I], from 
which we will construct all of the strong-segregation 
phases. Our calculations are performed in terms of the 
ratio 



a(A) 



_ A[z) 
~ A(R) 



(A = z/R) 



(2.2) 



of the cross-sectional area A at a height z relative to that 
of the outer surface, in an infinitesimal wedge of height 
R. This function may be easily calculated for wedges 
of particular shapes by elementary geometry, and is col- 
lected in Table H for various geometries. Since a(A) is the 
projected surface area along the normal vector extend- 
ing from the wedge point to the flat wedge top, it will 
be a quadratic function of A. The boundary condition 
a(0) = implies that a(A) is a sum of A and A 2 terms, 
and the boundary condition a(l) = 1 fixes the sum of 
the corresponding coefficients to be 1, leaving a single 
parameter. Hence we may generally write 



a(X)=p\+(l-p)X 2 



(2.3) 



The location Zd of the "dividing surface" separating 
the two species is determined by equating the relative 
volume below z, denoted v(z/R), to the volume fraction 
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f z/R Viz) 

<*/R)^l ^) = sm (2-4) 

v(p) = cf>v(l) (P = z d /R) , (2.5) 

where V(z) is the (partial) volume of the wedge below 
height z (see Figure |l|). 



TABLE I. Area function a(A) (see eq 2.2) for various wedge 
geometries, where A = z/R and R is the wedge height. The 
expression for p for the D, P, and G wedges refers to the 
geometry in Figure [l4|. 

Structure a (-M 
Lamellae 1 
Cylinders (p = 1) A 
Spheres (p = 0) A 2 
Symmetric Wedge (Figure 7]) (p = 2) A(2 — A) 

D, P, G Wedges (Figure 111) p A + (1 - p)\ 2 



p = h- (zi-z 2 ) x (d-b + a-c)/2n- (a-c) x (d-b) (D,P) 
n (zi-z 2 ) x (d-b)/2n- (a-c) x (d-b) (G) 
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FIG. 2. Area functions a(A) = 1, A — A 2 , A, and A 2 , for the 
lamellar, symmetric wedge, cylinder and sphere phases. 

The free energy per molecule of a wedge in the strong 
segregation limit is the sum of stretching and interfacial 
contributions, 



f fint ~t~ fstr ■ 



(2.6) 



The interfacial free energy per molecule is simply the sur- 
face tension contribution which, per chain, is 



fit 



7 Oo(g) 
v(l)R 



(2.7) 



where the surface tension 7 scales as y}l 2 '. j3^| In the 
strong segregation limit we ignore the translational en- 
tropy of the junction points, which scales logarithmically 
with molecular weight and is thus subdominant. 

The stretching energy is calculated by methods devel- 
oped for polymer brushes, ||^3)- The copolymer chains 
are added one by one, and the work to add each is 
summed. The height of the layer h A when the number of 
chains per area is a{h A ) is given by 

v(fi) - v(/3 - h A /R) = a{h A )n A n A a{p)/R (2.8) 
v((3 + h B /R) - v(f3) = a(h B )n B n B a(p)/R (2.9) 

where h A is the height of the growing A-layer, which is 
measured relative to the junction at z^', and similarly for 
B. 

An important quantity is the monomer chemical po- 
tential fi((f>) (the hydrostatic pressure, in an incompress- 
ible system), which is a decreasing function of distance 
from the dividing surface and is microscopically respon- 
sible for the stretching of the chains as their monomers 
seek regions of lower chemical potential. Under the as- 
sumption that there are free ends at all distances from 
the dividing surface, the chemical potential is quadratic 
in the distance from the dividing surface: p3| 



3^ 2 



(z A - z d y 



(2.10) 



A similar equation holds for the i?-species. For the in- 
wardly curved parts of the structure, this is exact; for 
the outwardly curved parts, this assumption leads to an 



unphysical negative density of free ends |J , but has been 
shown to give extremely good estimates of stretching free 
energy even for layers with curvature radii comparable to 
their thickness. |p4| 

The work to add an A-block is independent of the lo- 
cation of the free end, and so may be conveniently taken 
to be the work to add a chain with its conformation very 
near the surface, simply fl A ^i(zd). Hence the total free 
energy of a chain is obtained by integrating up to the 
desired coverage a, 



fstr — 



3tt^ 

8a 



da' 



■K (<? ) + 



Rl 



Rl 



h 2 B (a') 



Using eqs 2.8-2.9 we can write 



da = 
da = 



dh A a((3 + h A /R) 



dh f 



a((3) ' 
a(P-h B /R) 



fi(l - 4>) a((3) 



(2.11) 

(2.12) 
(2.13) 



Changing variables from a to h A and h B and using 
eqs 2.8-2~9|, we rewrite eq 2.11 



as 



fstr 



where 



In 



TT 2 R 2 

A =3 
3 



(1 - 4)R\ 



dya{(i-y)y 2 



1-/3 



dya(y + /3)y 



(2.14) 

(2.15) 
(2.16) 



At this point we make contact with previous calculations 
of asymmetric block copolymers and introduce an asym- 
metry parameter e: |2q,p2] 



K^Rl 

n\ Rl n A 



= —z-£ f 



(2.17) 



The second equality relates e 2 to Fredrickson's asymme- 
try parameter e F e.g. Ref. J3j|). The ratio t A = Q. A /R 2 A is 
a characteristic length which is independent of the length 
of an A-arm, and is larger for more flexible chains at a 
given volume. A smaller e indicates an enhanced ten- 
dency for the B species to stretch. Table [n] shows asym- 
metry parameters for several diblocks. 

TABLE II. Asymmetry para meter s for diblock copolymers 
(n A — n B — 1) as defined by eq [2. 17]. 



A-B 



ef = £ 



ref 



PE-PEP 
PE-PEE 
PEP-PEE 

PI-PS 



1.22 
1.58 
1.22 
1.27 
1.22 



1.5 
2.5 
1.5 
1.61 
1.5 



4 



1.11 1.23 |47] 

PB-PI 1.18 1.39 |47] 

PB-PS 1.31 1.72 47] 

PE=polyethylene, PS=l,4-polystyrene, PI=polyisoprene, 
PB=l,4-polybutadiene, PEP=att-poly(ethylenepropylene), 
PEE=poly(ethylethylene) 
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Upon minimizing eq 2.6 over the scale of the structure, 
i.e. the radius R, it is convenient to normalize all energies 
by a characteristic energy /o to obtain a compact form 
for the free energy: 



/ = (4 



where |52[ 



eh 



v{l) 3 \e(t) 2 (1 



1/3 



(2.18) 



s A = 



A a (R a )a a (/3 a ) 

a 

^v a {\)R a A a {R a ) 

OL 

^ R a A a {R a )I A (a, p a ) 
^v a {\)R a A a (R a ) 



(2.24) 



(2.25) 



1 /3 

Jo = (^) nV3 72 /3 (nAnB) l/3 ( ^ s) l/6 (219) 

is the free energy of a symmetric (e = 1) lamellar phase. 



B. Composite structures 

The procedure above applies to a single wedge. For 
the classical cylinder and sphere geometries with circu- 
lar (spherical) unit cells, each wedge is identical. For 
non-circular unit cells and for the complex geometries 
of bicontinuous phases, we must assemble the structure 
from many different wedges and minimize over the scale 
factor for the entire structure. The average free energy 
per chain f tot of a structure with many distinct wedges 
is 



fti 



E a UdV a 



V 



(2.20) 



where f a is the free energy per molecule in wedge a of 



volume dV a , and 



(2.21) 



is the volume of the structure. We choose a single scale 
factor i?o to determine the size of the whole structure. 
Each wedge a has its own area and volume functions 
a a (/3 a ) and v a (/3 a ), where the position of the dividing 
surface of each wedge, (3 a = Zda/Ra, is determined by 



V a (Pa) = <M*(1)- 



(2.22) 



The dimensionless functions {a a } are generalizations 
of eq 2.2 for each wedge a with wedge height R a . These 



functions encode the geometry of the particular struc- 
ture, and the cross-sectional area at the top of each 
wedge, A a {R a ), scales as Rq^ 1 , for a c?-dimensional struc- 
ture (e.g. d=l for lam ellae, d= 2 for cylinders, d = 3 for 
spheres). Expressing eq 2.7,2.14 in terms of {a ai A a } and 
Ro, we minimize over Rq to find the following free energy 
per chain of a particular structure: 



/ = /o S 2 a 



z<t> 2 (1 



1/3 



(2.23) 



where 



where I A (a,f3 a ) is obtained from eqs 2.15- 2.16] by sub- 
stituting a a (P a ) in place of a(/3), and a similar relation 
defines S B . 

By specifyi ng th e volume fraction in each wedge ac- 
cording to eq 2.22, we locally satisfy the constraint aris- 
ing from the fixed composition of the copolymers. In 
contrast, Likhtman and Semenov p2[| satisfied this con- 
straint only globally within a particular structure, which 
would be relevant for mixtures of different diblock copoly- 
mers with overall composition cf> |3(| in the strong seg- 
regation limit, in which the entropy of mixing of such 
different copolymers would be negligible 



III. 



FREE ENERGIES OF CLASSICAL DIBLOCK 
TOPOLOGIES 

A. Round unit cells 



Using the results of Section II A we can find the ener- 
gies of the classical phases of diblock-copolymers: lamel- 
lae (L), cylinders (C), and spheres (S) in the round unit 
cell approximation, in which the unit cells are taken to 
consist of identical wedges. The corresponding free ener- 
gies are: 



flam 
fey I 
fsph 



el 



1/3 



2^(1-^/2)3(3 + 01/2) 20 

(1 - <t>) 2 + T 



-,1/3 



= 3 



•0 4 / 3 (l-0 1 /3)3^2/3+3^1/3 + 6 ) ^ 



10(1 -<pf 



10e 



1/3 



(3.1) 
(3.2) 

3 

(3.3) 



Calculations based on round unit cells || provide lower 
bounds for the free energy, because they in fact describe 
the free energy per molecule of micelles. M We may imag- 
ine a volume packed with such micelles, the interstitial 
regions filled with compatible long homopolymer with 
negligible surface tension against the outsides of the mi- 
celles, and negligible entropy of mixing. Then we could 
do work to deform the micelles into a space-filling array, 
expelling the homopolymer at no free energy cost. To 
distinguish between crystal structures within a particu- 
lar topology, such as between hexagonal and square for 
the cylindrical topologies, we must examine the energy 
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for packing the molecules into the particular geometry, 
which is performed below. 



1 .2 




FIG. 3. (Above) Hexagonal and round unit cells for the 
cylindrical phase, and (below) typical wedges for spherical, 
cylindrical, and bicontinuous phases; the dividing surface is 
shaded. 
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FIG. 4. Free energies of classical phases and the symmet- 
ric wedge, for straight-paths constructions (asymmetry factor 
e — 1). The shaded regions are bounded by the upper (BCC 
and HEX, respectively) and lower (spherical and cylindrical 
micelles, respectively) free energy bounds for the spherical 
and cylindrical topologies. The symmetric wedg e has area 
function a(A) = 2(A — A 2 ), described in Sec. [V A, 



B. Non-round unit cells (straight paths) 



We can produce an upper bound for different struc- 
tures by assembling small pieces of the cylindrical or 
spherical micelles to fill the appropriate unit cell. Each 
wedg e has a parabolic monomer chemical potential given 
by eq 2.10 . However, each wedge a has a slightly different 
shape and geometry, and thus has a distinct potential fx a . 
Adjacent wedges are not in equilibrium with each other 
and will relax if allowed to do so. Hence the calculation 
yields an upper bound. To construct the unit cell of, e.g., 
hexagonal cylinders, we assume a hexagonal dividing sur- 
face scaled down by cj) 1 / 2 and assemble the unit cell from 
tiny pie-shaped wedges extending from the center of the 
hexagon to the cell boundary. We make an analogous 
construction for square arrays of cylinders, or for FCC 
and BCC packings of spheres. 

We calculate the volume- avera ged stretching free en- 
ergy per molecule using eq 2.23 . To calculate this in 
practice we use the following procedure. For cylindrical 
micelles we divide a cell of a given symmetry (say, hexag- 
onal) into tiny wedges. Each wedge is adjusted slightly 
by making the segment of the wedge on the dividing sur- 
face normal to the bisector of the wedge, which is the 
path of the polymer. Such an adjustment introduces a 
negligible volume in the continuum limit of many small 
wedges. T he sur face area used for calculating the surface 
energy (eq 2.24) is, of course, the area of the segment in 
the original hexagonal dividing surface (before adjusting 
the wedge to account for straight paths). 



For the classical phases (see Figure ^) the ratios of the 
upper and lower bounds for the free energies are inde- 
pendent of <j), given by: 



fhe 



fey I 
fsquc 



10 



cyl 



1/3 



1/3 



fb, 



fsph 

[fee 
fsph 



95 \ 1/3 fl 



384 J V 2 



V3 



2/3 



1/3 



1.036 



1.101 



1.072 



1.077. 



(3.4) 



Evidently, the most favorable structures have the 
"roundest" unit cells. The hexagonal phase is favored 
over the square phase, and BCC is slightly favored over 
FCC. 
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C. Non- round unit cells (kinked paths) 




FIG. 5. Hexagonal unit volume divided into wedges for 
straight-path (left) and kinked-path (right) ansatzes, for equal 
compositions (j>. Also indicated are adjustments necessary to 
calculate the free energy in each wedge for the straight-path 
ansatz. 



and similarly for dR' . Remarkably, the upper bound for 
the kinked-path-hexagonal ansatz is typically less than 
1% above the lower bound of cylindrical micelles, and the 
transition is shifted to only a slightly smaller A fraction 
</> (see Figure ^). Apparently the extra stretching energy 
to maintain a hexagonal A — B interface with straight 
paths is relaxed considerably by allowing the inner block 
to adopt a more nearly circular dividing surface, which is 
preferred. Recent accurate numerical self-consistent field 
calculations of diblock melts have shown that in fact the 
A— B interface is nearly circular, with a slight hexagonal 
modulation (angular modulation with 6-fold symmetry) 
of relative amplitude 0.03% at xN = 60 and <fi — 0.33 




Other upper bounds can be obtained by using differ- 
ent prescriptions for the A—B surface. For example, one 
could choose a circular A — B surface of radius D(<p) for 
the cylindrical phase. The inner (A) volume may be di- 
vided into wedges, and the outer (B) volume divided into 
wedges which each satisfy the volume constraint <f> with 
a partner A wedge |2(|. For a right triangular wedge 
which subtends an angle 6*o, points at angle 9 on the 
A—B surface map to points R(0) on the boundary of the 
Wigner-Seitz cell, 



R(0) =Ri + a(0) (R2-R1) 



(3.5) 



where Ri = {0, 1}, and R2 = {tan#o, 1}. The composi- 
tion specifies the radius, according to 



> tan O = D' 



(3.6) 



The mapping which obeys the local composition con- 
straint is 



<0) = 



D(<j>) 



sin 9 
tan 6*o 



1 - D{4>) cos( 



(3.7) 



where 9o = tt/3 and 7r/4 for hexagons and squares (see 
Figure ||) , respectively. 

Minimizing over the scale of the structure, eq 2.20] 
yields, after some calculation, the following free energy: 



f_ 

fo 



where 



2c6 



(1 - ^D^<f) 



d9 



R 



dR' 




dr' 






d9 


+ 


d9 


i 


) 



1/3 



r{9) = D(4>) 
r' = r ■ 6 



sin 9 
cos 9 
(R-r)(R-r) 



IR 



(3.8) 
(3.9) 

(3.10) 
(3.11) 
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FIG. 6. Free energies for cylindrical topologies in var- 
ious approximations (asymmetry factor (e = 1). Also 
shown is the free energy of the symmetric wedge. The in- 
set shows an enlarged view of the crossing of the wedge 
and kinked-path-hexagonal free energies. The lamellar en- 
ergy crosses the cylindrical micelle energy at (f> = 0.299, 
the kinked-path-hexagonal energy at cj> = 0.293, and the 
straight-path-hexagonal energy at cj> — 0.255. 



IV. GEOMETRY OF BICONTINUOUS PHASES 



A. Generic saddle surfaces 



Before addressing particular symmetries (P, D, or G) 
of bicontinuous phases, we discuss the closest analogue to 
a round unit cell. We would like to produce a simple esti- 
mate of the free energy, analogous to the cylindrical and 
spherical micelle calculation, which captures the physics 
of bicontinuous topologies. We thus represent a generic 
bicontinuous phase as a wedge, shown in Figure |7]: an 
infinitesimal patch of "saddle" surface, with edges given 
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by the normals, terminating in a small line segment ly- 
ing along the bond-lattice. We envision the surface as a 
minimal surface, which has zero mean curvature. pq] 



The stretching free energy per molecule of the symmet- 
ric wedge can be calculated as before, given the relative 
area as a function of relative height along the center nor- 
mal (Table @): a(A) = A(2 - A), where A = z/R, z = is 
the thin end of the wedge and z = R is the patch of mini- 
mal surface. As before, the dividing surface location Zd is 
determined by eq 2.5. The resulting free energy is given 
by applying eqs |2.5| -2.7, and is shown with the various 
cylindrical bounds in Figure |^. This estimate misses by 
a few tenths of one percent the intersection of the lamel- 
lar phase and the kinked-path upper bound bound for 
the hexagonal phase, and is stable with respect to the 
straight-path upper bound. 




z=z. 



z=0 



FIG. 7. Symmetric saddle chip (p = 2 in Eq.gJ). The A—B 
dividing surface at z = Zd is shaded. A material fills the top 
and bottom of the wedge, and B material fills the volumes 
near the (minimal) partitioning surface. 



Clearly the simple wedge construction captures some 
important physics. The structures formed by copolymers 
at different volume fractions <f> arise from competition 
between interfacial and stretching free energies. The dif- 
ferent structures present different functions a(A), which 
determine both the dividing surface area and the stretch- 
ing energy as a function of volume fraction. The phases 
occur in the order they do because the progression of 
functions a(A) from quadratic A 2 (spheres) to linear A 
(cylinders) to A(2— A) (bicontinuous) to constant (lamel- 
lae) gives progressively less volume to the "outer" chain 
to avoid stretching, but uses progressively less area to 
separate the two species at higher volume fractions of 
the minority species. 
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FIG. 8. Free energies for saddle wedges of various shape 
factors p — 1.5, 1.9, 2.0, 2.1. Also shown is the free energy of 
the kinked-path hexagonal phase. 

However, we cannot argue as before that our simple 
estimate for the bicontinuous phase is a lower bound for 
actual bicontinuous phases, because there is no way to 
pack together copies of one infinitesimal wedge to pro- 
duce a "micelle" that 1) fills some region of space, and 
2) is bounded by some surface(s), the volume outside of 
which could be filled by homopolymer. Neither can we 
argue that this estimate is an upper bound, because we 
certainly cannot pack a unit cell of the region bounded 
by the D or G surfaces with identical copies of one in- 
finitesimal wedge. 

Bicontinuous phases are assembled from different 
wedges a, with different shape factors p a in Table Q. 
The shape factor roughly gauges the splay or Gaussian 
curvature of the surface at the top of the wedge, with 
p = 1 (cylinders) corresponding to zero Gaussian curva- 
ture, p < 1 to positive Gaussian curvature, and p > 1 
to negative Gaussian curvature. [The Gaussian curva- 
ture is the product of the two radii of curvature of a 
surface] . The distribution of wedges must be chosen to 
pack the desired structure. While the symmetric wedge 
has p = 2, this is not an optimum shape. In fact, the 
optimum wedge shape depends on composition, as can 
be seen in Figure |[ It is evident that there are shape 
factors p which have lower free energy than the straight- 
and kinked-path hexagonal upper bounds (Figure ^), so 
it is not unreasonable to hope that a judicious packing 
configuration can be a stable thermodynamic phase. 

The effect of the shape factor p on the topology of the 
phase diagram emerges upon examining conformation- 
ally asymmetric (e ^ 1) copolymers. Following Ref. p5[ , 
we explore the effect of conformational asymmetry on 
the stability of bicontinuous phases by multiplying the 
wedge free energy by an additional arbitrary small pref- 
actor (0.99) which enhances stability. [We will see be- 
low that for e = 1 the bicontinuous phases that we can 
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calculate (G, D, P) are of order 0.3 — 2.2% higher than 
the cylinder-lamellar crossing, depending on which upper 
bound one compares.] Figure [To| shows 'phase diagrams' 
as a function of conformational asymmetry e. 

Recall that for e > 1 the _B-block is more flexible, 
while the A-block is stiffer and better able to stretch. 
For the symmetric wedge (p = 2) conformational asym- 
metry reduces the stability of the stiff-minority wedge 
phase (4> < 1/2) and enhances the stability of the flexible- 
minority wedge phase (4> > 1/2), and shifts all transitions 
to greater <f). For p > 2.0 the wedge phases lose stability, 
as could be guessed from Figure |[ For p < 2.0, for which 
the wedge is more cylindrical-like (p = 1 coresponds 
to cylinders), conformational asymmetry enhances the 
stiff-minority wedge phase relative to both the lamellar 
and cylindrical phases, and decreases the stability of the 
flexible-minority wedge phase. We emphasize that these 



p=1.8 








are not phase diagrams, for a true phase is a mixture of 
wedges with different shape factors which fill space. 
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FIG. 9. Free energy relative to that of the kinked-path 
hexagonal phase, as a function of saddle wedge shape factor 
p, for compositions (j> — 0.28, 0.31, 0.35. 



p=1.9 




FIG. 10. Phase diagrams for the round unit cell approximation, including the saddle wedge for shape factors 
p — 1.8, 1.9, 2.0, 2.1, as a function of composition (j> and conformational asymmetry e. The free energy of the wedge has been 
multiplied by an additional prefactor (0.99) to explore the effect of conformational asymmetry on the stability of bicontirmous 
phases. 

B. Conformational Asymmetry 

Experimentally, the G phase has been observed be- 
tween the lamellar and cylindrical phases in several 
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strongly-segregated copolymer systems, for <fi around 0.3 
(and, symmetrically, around 0.7). Some groups have ar- 
gued for phase stability of bicontinuous phases in terms 
of bending rigidities, |p9| , ^o|| . However, there is funda- 
mentally no bending energy in the problem; descriptions 
in terms of bending energies only arise from a proper 
accounting of stretching free energy in curved geome- 
tries. Our approach is to choose a geometry as an ansatz 
and compute the corresponding interfacial and stretch- 
ing free energies in a manner consistent with calcula- 
tions for cylindrical, spherical, and lamellar phases. The 
structure, revealed by scattering and electron microscopy, 
pl| , ^3[ studies, can be described as follows, p8|Jtl| |. 



Consider first the D geometry, which is easier to vi- 
sualize. A skeleton formed of the bonds of a diamond 
lattice is shown in Figure O. Two such lattices interpen- 
etrate, analogous to the interpenetration of two simple 
cubic lattices in a BCC structure. Now imagine swelling 
the bonds in these lattices into tubes of a finite diame- 
ter. The walls of these tubes are a rough approximation 
of the experimentally observed "dividing surface" sepa- 
rating the regions containing the two blocks. The volume 
contained within the tubes corresponds to the region in- 
habited by the low volume-fraction monomer. 



To model the D geometry, we use a self-dual minimal 
surface, called the Schwartz D (diamond) surface, which 
partitions space into two identical interpenetrating re- 
gions, each of which contains and is topologically equiv- 
alent to a diamond bond-lattice ^|l(],[ll]] . Within each 
of these regions is a dividing surface, which surrounds a 
copy of the bond- lattice. The copolymer chains then have 
conformations with one (A) species stretching towards 
the bond-lattice, the junction between blocks residing on 
the dividing surface, and the other (B) species stretch- 
ing towards the minimal surface. In the G phase the 
diamond lattice is replaced by a three-fold coordinated 
lattice, and the surface is replaced by the gyroid minimal 
surface discovered by Schoen in 1970 fl4[ ], in which the 
two interpenetrating G volumes are chiral enantiomers of 
one another. For the P phase the bond lattice is six-fold 
coordinated (Figure 11) and the candidate partitioning 
surface is the Schwartz P minimal surface. 



There is no compelling reason to choose a minimal sur- 
face for the partitioning surface. However, minimal sur- 
faces solve the variational problem of minimizing surface 
area with zero pressure across the interface. If the diblock 
phase is in fact partitioned into two equivalent connected 
regions, then by symmetry there can be no net pressure 
exerted across the dividing surface that separates the two 
equivalent disjoint connected regions. So a minimal sur- 
face is reasonable, but by no means certain, since there 
is no obvious area energy to minimize. 




FIG. 11. Skeleton bond lattices for, from top to bottom, 
the P, D, and G surfaces. 



The P, D, and G surfaces may be conveniently calcu- 
lated using the Weierstrass representation. jj8| Here, the 
three-dimensional points r of the two-dimensional surface 
are parametrized by the complex number to — u>i + zaj 2 . 
The G, P, and D surfaces are triply-periodic minimal 
surfaces with space groups laid, Im3m, and Pn3m, ve- 



il 



spectively. To generate the full surface it is enough 
to calculate a single patch, to which all the symmetry 
operations of the space group may be applied to gener- 
ate the full structure. A generic surface has two radii 
of curvature which are generally non-zero and different. 
Points where the surface is flat are singular points, since 
both radii of curvature are zero and a direction of the 
surface cannot be determined. These flat points define 
the corners of the fundamental patch. 
The Weierstrass representation is: 



Re 



v(w) e 



(w 8 - 14cj 4 + 1)V2 



where 



v(«) - {1 



u 2 ,i i 



1 



(4.1) 



(4.2) 



and V is the domain of integration shown in Figure |l2j. 
The points on the corners correspond to the flat points, 
and it is evident that the integrand above (excluding the 
measure) is singular at these points. The angle 9 deter- 
mines the surface: 




D (Pn3m) 
P (7ro3m) 
G {laid) 



(4.3) 




FIG. 12. The thick line bounds the domain of integration 
V for the Weierstrass representation. The arcs are from circles 
of radius y/2 centered at the points (±1, ±l)/v2. 



FIG. 13. Top: fundamental patches of, from left to right, the P, D, and G surfaces. Below: larger portions of the surfaces. 



For other angles the surface intersects itself. This does 
not, of course, exhaust the class of triply periodic min- 



imal surfaces either mathematically, |14[|41|jl2|| or phys- 
ically 43 1 . We have chosen the D, P, and G surfaces 
because they are the most common observed surfactant 
bicontinuous surfaces, and have been claimed experimen- 
tally in block copolymers. Sections of these surfaces are 
shown in Figure 13 For fairly accurate calculations (yield- 
ing energies lower than those for the true surface by of 
order a few tenths of a percent) the D surface may be ap- 



proximated by a simple hyperbolic surface, z = xy. This 
suggests that the minimal D surface may not the opti- 
mal partitioning surface. However, we have varied the 
shape of the partitioning surface around the D surface, 
and found free energy variations of only a few tenths of a 
percent. Because of this, we have not optimized the free 
energy with respect to adjustments in the partitioning 
surfaces. 

To produce an upper bound on the free energy of 
the bicontinuous phases we follow a procedure analogous 
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to that used for the classical cylindrical and spherical 
topologies. Namely, we divide a unit cell of these struc- 
tures into a large number of wedges (shown in Figure 14), 
similar to Figure 1, but of varying radii and Gaussian 
curvat ure, and a verage the free energy per molecule us- 
ing eq 2.23 2.25] . Each structure has a different unit cell 
(a big wedge) from which the entire structure may be 
generated by applying the symmetry operations of the 
particular space group. Figure [l5] shows the fundamen- 
tal cell for the G structure. 

This procedure is straightforward. First we calculate 
the minimal surface, and then adopt a convenient map- 
ping from the points on this surface to the skeleton. Es- 
sentially, we construct an interpolation between those 
high-symmetry points on the minimal surface with nor- 
mals that project onto the underlying bond-lattice. For 
the D and G phases, we optimize the mapping from the 
partitioning surface to the line segments of the bond lat- 
tice to minimize the free energy. We perform this by 
a conjugate gradient algorithm that distorts the two di- 
mensional mesh of points on the surface, and gains of 
order 1% in energy. 

ft 




FIG. 14. Fundamental wedges used for constructing the D 
or P (left) and G (right) structures. We have shown quite gen- 
eral wedges, before they undergo slight adjustments to ensure 
that n is normal to the bond segment zi — Z2 and the surface 
element abed. 



line segment connecting the center of the outer surface 
of the wedge to the center of the hexagon. Such adjust- 
ments are negligible in the limit of infinitesimal wedges. 
As before, the area of the A— B interface is the true area, 
rather than the (smaller) area that results from adjust- 
ing the wedge to assure a chain path normal to the A— B 
interface. 

In this way, the unit cell of the region bounded by a 
minimal surface is decomposed into many small wedges, 
each with a known (and different) radius R, shape fac- 
tor p, and volume. The locatio n of the dividing surface 
within each wedge is fixed by eq 2.22| , and the free energy 
calculated with eqs 2.23-2.25. We have checked the algo- 



rithm and the dependence on the fineness of the mesh by 
using it to successfully compute the free energies of the 
hexagonal phase. 



V. RESULTS 
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FIG. 15. Side and top views of basic unit for G structure, 
showing the initial division into wedges. 

Thus, each small patch on the minimal surface is con- 
nected by straight lines to a small line segment on the 
skeleton, and a set of wedges results. Each wedge is ad- 
justed slightly by making both the top patch and the 
bottom segment orthogonal to the line segment connect- 
ing the center of the patch to the skeleton. This is anal- 
ogous to making the outer surface of the wedge in the 
upper bound for the hexagon phase orthogonal to the 



FIG. 16. Free energies for bicontinuous phases (upper 
bounds) compared to upper and lower (cylinders) bounds of 
the hexagonal phase. The inset enlarges the crossing of the 
G and straight-path-hexagon upper bounds. Asymmetry pa- 
rameter e = 1.0. 

Figure [l^ shows the free energy curves as a function 
of <p for conformational^ symmetric copolymers (e = 1). 
The upper bound for the ImZm phase (P) lies 4 — 5% 
above that for the Pn3m phase (D), which in turn is 
less than a percent (~ 0.7% at (f> = 0.3) above the 
IaZd (G) phase. Consistent with experiments and self- 
consistent field theory, we do not find a stable G phase. 
At the lamellar-kinked-path-hexagons free energy cross- 
ing {(p ~ 0.293) the free energy of the G phase is of 
order 2.2% larger, while at the lamellar-straight-path- 
hexagons free energy crossing (</> ~ 0.255) the free en- 
ergy of the G phase is only a few tenths of a percent 
(~ 0.28%) greater. This, we have argued, may be the 
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fairer comparison, since both calculations use straight 
paths. Unfortunately, we do not know how to perform a 
kinked-paths estimate for the bicontinuous phases. 

Note, however, that we do not expect to gain as much 
energy from a kinked-path calculation for the bicontinu- 
ous phases as for the cylindrical topologies. Consider the 
hexagonal calculation. The boundaries of the Wigner- 
Seitz cell, and hence the A — B dividing surface, have 
sharp corners into which the chains must stretch. Pre- 
sumably a large part of the gain in the kinked-path calcu- 
lation comes from relieving the strain associated with this 
stretch, and relaxing the inner block to its preferred cir- 
cular structure. Bicontinuous phases, on the other hand, 
have smooth "Wigner-Seitz boundaries" (i.e. the mini- 
mal surface), and expensive stretching occurs mainly at 
the junctions of the skeleton lattice. Hence, the anoma- 



lous stretching that may be relieved by a kinked-path 
calculation occurs along points in the structure, rather 
than along lines. So we expect that our straight-path es- 
timate is not likely to differ greatly from a kinked-path 
estimate, and that the free energy of the G phase re- 
mains well above that of the kinked-path hexagons. The 
stretching of the chains at the junctions presumably con- 
tributes to the relative stability of the P, D, and G struc- 
tures, which have 6-, 4-, and 3-fold coordinated bond lat- 
tices. The more highly-coordinated lattices require more 
chain stretching to accommodate the space, which sug- 
gests that P, D, and G occur in increasing order of sta- 
bility. This picture is corroborated by recent work of 
Matsen and Bates Q] , who quantitatively examined the 
packing frustration in the G, D, and HPL phases. 




FIG. 17. Free energy crossings for gyroid, hexagonal (kinked paths) and lamellar phases for asymmetry parameter values 
e = 0.7, 1.0, and 1.3. For e > 1 the species (B) with composition 1 — is more flexible and/or has more arms than the other 
(A) species, and the A species has an enhanced tendency to stretch. 



Our results apply to the strong-segregation limit, 
which is attained in the limit of large x- Because the 
phase boundaries shift away from (f> — 0.5 as x increases 
from weak segregation we expect our phase bound- 
aries to be further from (j) = 0.5 than experimental values, 
which is indeed the case. 



Figure [L7| shows free energy crossings for various values 
of the conformational asymmetry parameter e. Relative 
to a conformationally symmetric melt, conformational 
asymmetry stabilizes phases with a stiff minority species 
and destabilizes phases with a stiff majority species, mov- 
ing boundaries to larger 4>(A) for e > 1. Recall that for 
e < 1 the inner A-block is more flexible, while the B- 
block is stiffer and better able to stretch. We find that 
e < 1 reduces the relative stability of the G phase, while 
the stability is enhanced for e > 1, and becomes stable 
for rather large asymmetries e > 9.0 (Figure |I§|). 
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FIG. 18. Free energy of G phase relative to that of the 
kinked-path upper bound estimate of the lamellar-hexagonal 
phase boundary, as a function of conformational asymmetry 
e. The vertical line marks the conformationally symmetric 
copolymer, and the horizontal line is for reference. 

Previous calculations of phase diagrams of 
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conformationally-asymmetric diblocks have been done in 
the weak segregation regime, [fl5]]4(| and, for the generic 
symmetric wedge as a model bicontinuous structure, in 
the strong segregation regime. [E5| Matsen and Bates 
p6| found, as we do, that conformational asymmetry 
stabilizes the G phase with a stiff minority phase, widen- 
ing the composition window and moving the lamellar-G 
and G-cylinder boundaries to higher stiff compositions; 
and destabilizes the G phase with a stiff majority phase, 
both narrowing the composition window and shifting it 
to a higher stiff fraction. Their calculations are limited 
to xJV < 30, and it is inconclusive whether the %iV — * oo 
limit of this calculation yields a stable G phase. The 
same qualitative behavior was found for the generic sym- 
metric wedge in the strong segregation regime. [^5| 

Table [5] summarizes asymmetry parameters from re- 
cently collected data J47],fii|. While none of these di- 
blocks have the large conformational asymmetry required 
to test our prediction of a stable phase in strong segrega- 
tion, A n — B m starblock copolymers have asymmetry fac- 
tors larger, by a factor of n/m, than those due to intrin- 
sic chain stiffness effects alone. For example a PEg-PEEi 
starpolymer has an asymmetry factor e ~ 6x 1.58 = 9.48. 

We have attempted to calculate energies for the HPL 
p8| , p9]| phase, which is now accepted as a metastable 
phase piJI . The minimal crystal phases D, P, and G 
have obvious candidate minimal surfaces to act as an in- 
termatcrial dividing surface towards which the majority- 
phase ends stretch; and the minority-phase ends stretch 
towards the skeletal bond lattice. On the other hand, the 
majority-phase ends in the HPL phase stretch towards a 
combination of lines (in the hexagonally-arranged per- 
forating tubes) and surfaces (within the majority-phase 
layer). Similarly, it is not obvious how to partition the 
minority-phase ends between lines and surfaces. The re- 
sult is a non-analytic mapping which is difficult to mini- 
mize over. Our attempts have thus far yielded quite high 
energies, of order that of the P phase. 



is locally equal to <p. The interfacial contribution 
to the free energy is the area of the A-B interface 
times the A-B surface tension. 

4. Compute the stretching free energy per chain for 
each wedge within the approximation of straight 
paths. Calculations for straight paths involve slight 
adjustments to the wedges whose contributions 
vanish in the limit of small wedges. 

5. Optimize the free energy per chain with respect to 
the overall scale of the mesophase (e.g., the dimen- 
sion of the unit cell). 

6. Optimize the mapping from the partitioning sur- 
face to the bond skeleton to minimize the overall 
free energy. 

In certain structures (e.g. HPL) the majority phase ends 
lie on both lines and surfaces, in which case the procedure 
above must be suitably generalized. We have also shown 
how to calculate the free energy for geometries where the 
shape of the A — B interface is specified, for phases of 
cylindrical topology. This requires polymer chain paths 
which are kinked at the A— B interface. 

The infinitesimal wedg es a re described by the relative 
area function a(z/ R) ( eq 2.2) which is parametrized by a 
single scalar p (eq 2.3 and Table |) that roughly gauges 
the local Gaussian curvature of the partitioning surface. 
For the classical phases all wedges are identical, while 
bicontinuous phases have different distributions of shape 
factors p. For symmetric stars we find a metastable bi- 
continuous (gyroid, or G) phase which is most stable 
near the lamellae-hexagonal cylinder transition. For suf- 
ficiently asymmetric copolymers (e > 9.0) we predict a 
stable G phase. 



VI. SUMMARY 

We have outlined a general method for computing the 
free energy of block copolymer phases in the strong seg- 
regation regime. The procedure consists of the following 
steps: 

1. Choose a candidate geometry and an associated 
partitioning surface that divides space into disjoint 
interpenetrating regions (the majority blocks from 
the two regions stretch towards this surface). 

2. Divide the enclosed volume into infinitesimal 
wedges, defined by straight paths connecting the 
partitioning surface to a skeleton of bonds (the mi- 
nority blocks stretch towards this bond skeleton). 

3. The A-B interface in each wedge is located such 
that the fraction of wedge volume filled by A blocks 
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